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Abstract 

We propose a general framework for the simultaneous modeling of equity, govern- 
ment bonds, corporate bonds and derivatives. Uncertainty is generated by a general 
affine Markov process. The setting allows for stochastic volatility, jumps, the possi- 
bility of default and correlation between different assets. We show how to calculate 
discounted complex moments by solving a coupled system of generalized Riccati equa- 
tions. This yields an efficient method to compute prices of power payoffs. European 
calls and puts as well as binaries and asset-or-nothing options can be priced with the 
fast Fourier transform methods of Carr and Madan (1999) and Lee (2005). Other Eu- 
ropean payoffs can be approximated with a linear combination of government bonds, 
power payoffs and vanilla options. We show the results to be superior to using only 
government bonds and power payoffs or government bonds and vanilla options. We 
also give conditions for European continent claims in our framework to be replica- 
ble if enough financial instruments are liquidly tradable and study dynamic hedging 
strategies. As an example we discuss a Heston-type stochastic volatility model with 
possibility of default and stochastic interest rates. 

Key words Pricing, hedging, affine models, stochastic volatility, jumps, default. 

1 Introduction 



The goal of this paper is to provide a flexible class of models for the consistent pricing 
and hedging of equity options, corporate bonds and government bonds. The noise in our 
models is driven by an underlying Markov process that can generate stochastic volatility, 
jumps and default. For the sake of tractability we assume it to be affine. Then discounted 
complex moments of the underlying can be calculated by solving a coupled system of 
generalized Riccati equations. This yields an efficient method to compute prices of pov^^er 
payoffs and the discounted characteristic function of the log stock price. From there, 
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prices of vanilla options as well as binaries and asset-or-nothing options can be obtained 
with fast Fourier transform methods a la Carr and Madan |2] and Lee 1 15]. We also give 
conditions for European contingent claims in our models to be replicable if there exist 
enough liquid securities that can be used as hedging instruments. 

Our framework can be seen as an extension of the unified pricing and hedging model 
of Carr and Schoutens [3], where market completeness is achieved through continuous 
trading in the money market, stock shares, variance swaps and credit default swaps. 
The authors suggest to approximate general payoffs with polynomials. In this paper we 
propose an approximations with a linear combination of government bonds, non-integer 
power payoffs and European calls. 

Affine models have become popular in the finance literature because they offer a good 
trade-off between generality and tractability. One-factor affine processes were first used 
by Vasicek Il22ll and Cox-Ingersoll-Ross Hi for interest rates modeling. Popular affine 
stochastic volatility models include the ones by Stein and Stein (2T\ and Heston [9]. For 
affine models in credit risk we refer to Lando [|13.I . Here we work with general affine 
processes in the sense of Duffie et al. O. 

The rest of the paper is organized as follows. In Section 2 we introduce the model. In 
Section 3 we show how discounted complex moments can be calculated by solving gener- 
alized Riccati equations. This yields an efficient way of calculating options with power 
payoffs. As a corollary one obtains conditions for the discounted stock price to be a mar- 
tingale under the pricing measure. The prices of vanilla options as well as binaries and 
asset-or-nothing options can be computed with the fast Fourier transform methods of Carr 
and Madan [2] and Lee [15 |. For the pricing of European options with general payoffs we 
propose an -approximation with a linear combination of government bonds, power pay- 
offs and European calls. We illustrate this method by pricing a truncated log payoff, which 
can be applied to the valuation of a variance swap in the case where the underlying can 
default. Section 4 is devoted to the derivation of hedging rules. We show that in a model 
with discrete jumps, or no jumps at all, every European option can perfectly be hedged 
by trading in stock shares and a proper mix of European calls, government and corporate 
bonds. A system of linear equations is derived to find the hedging strategies. As an exam- 
ple we discuss in Section 5 a Heston-type stochastic volatility model with the possibility 
of default and stochastic interest rates. 

2 The model 

Let {Xt,Fx)t>o, xeD be a time-homogeneous Markov process with values in D := R™ x 
M"""*, m G N, n G N \ {0}. Let be the expectation corresponding to and denote by 
(•,-) the Euclidean (non-Hermitian) scalar product on C": {x,y) := Y17=i^iyi' ^^V ^ C*^. 
Furthermore, set Z = {1, . . . , m} and J = {m + 1, . . . , n}. Xx^t and Xj^t denote the first m 
and last n — m components of Xt, respectively. 

We model the risk-neutral evolution of the price of a stock share by 

St = exp(st + Rt + At)l{i<^}, 
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where 



St = e + (e, Xt) for (e, e) G M x M" 

Rt = f Tudu for rt = d+ {6, Xx,t) , {d, S) € I 
Jo 

At = Kdu for At = c + (7, Xj,t) , (0,7) G 

JQ 



and 



r = inf{t > : At = ^} 



for a standard exponential random variable E independent of X. 

The process rt models the instantaneous risk-free interest rate and At the default rate. 
Note that both of them are non-negative, r is the default time {St = for t > r), and st 
describes the excess log-return of St over rt + At before default. Alternatively, one could 
model the three processes st = st+Rt+At, Rt and At. But modeling (st, Rt, At) is convenient 
since it makes the discounted stock price equal to exp(st + At)l{t<r} and will allow us to 
give a simple condition for it to be a martingale under P^; (see Corollarv 13.51 below). 

For the sake of tractability we always assume that {Xt,^x)t>o,x(^D is stochastically 
continuous and affine in the sense that 

(HI) Xt — > in probability for t ^ to with respect to all P^;, x e D, and 
(H2) there exist functions : M+ x iP" ^ C and V : P+ x zP" ^ C" such that 

Kx[ex.p{{u, Xt))] = exp + {ip{t,u),x)) 

for all X G -D and {t, n) G P+ x iP". 

It has been shown by Keller-Ressel et al. lfT2]| that under the hypotheses (HI) and (H2), 
{Xt,fx)t>o, xeD automatically satisfies the regularity condition of Definition 2.5 in Duffie 
et al. O. Therefore, according to Theorem 2.7 in Duffie et al. [6], it is a Feller process 
whose infinitesimal generator A has C^{D) (the set of twice continuously differentiable 
functions f : D with compact support) as a core and acts on / G C^iD) like 



+ {b + /3x,S/f{x)) 



k,l=l 
+ 



Oxkdxi 

{fix + - fix) - >yjfix),xjim Hdo 



(2.1) 



D\{0} 
m „ 

■ E / (/(^ + ^) - /(^) - {^JU{^}fix),XJU{^}iO)) X^^^iidO, 

Jd\{o} 



where x = (xi,---,Xn) for 



XkiO 



.(lAlal)fr 



otherwise 



and 

(i) a G 



is symmetric and positive semi-definite with axx = 
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(ii) a = {ai)i^x, where for each i £ I, at £ M"^" is a positive semi-definite symmetric 
matrix such that all entries in Oj xx are zero except for ai^u 

(iii) beD 

(iv) /3 G M"^" such that (3xj = and j3xx has non-negative off-diagonal elements 

(v) is a Borel measure on -C>\{0} with 

/ ((xx(o,i) + iix^(e)in^K)<oo 

Jd\{o} 

(vi) fi = {fii)i£x is a vector of Borel measures on I?\{0} such that 

/ {{XX\{^}{0, 1) + \\XJU{^}m\^) l^MO < OC. 

Jd\{o} 

In addition we require the jumps to satisfy the following exponential integrability condi- 
tion: For all q G W'\ 

(vii) 

/" e<5'«h|||^ll>i}z/(d0 < and /" e<«'«h|||^||>i}^i(de) < oo for alH G X. 
Jd Jd 

It follows from Theorem 2.7 and Lemma 9.2 in Duffie et al. [6] that for every set of 
parameters (a, a, 6, /3, ji) satisfying (i)-(vii), ( 12. ID defines the infinitesimal generator of 
a D-valued Feller process iXt,Fx)t>o, xgd satisfying conditions (HI) and (H2) for {(pjip) 
equal to the unique solution of the following system of generalized Riccati equations 

'dt(l){t,u) = Fo{i;{t,u)), 0(O,u)=O 

< dtMt,u) = F{i;{t,u)), M'^,u) = ux , (2.2) 

where the functions Fq : C"- ^ C and F : C" ^ C™ are given by 

Fo{n) = {an, u) + (6, u) + f (e<"'«> - 1 - {uj^xAO)) HdO 

Jd\{o} ^ ' 

n 

Fi{u) = {aiU,u) + \2l^kiUk+ [e^'''^^ -I- {uju{i},Xju{i}{0)) fJ-iidO, i^'^- 

In particular, (a, a, 6, /?, fi) uniquely determine the transition probabilities of X and 
therefore, its distribution under each ¥x, x £ D. Since every Feller process has a RCLL 
version, we may assume (X* , r, P:r)t>o,xe d to be defined on TZd x M+, where TZd is the space 
of all RCLL functions uj : M+ — > D, and (Xt, r){uj, y) = {uj{t),y). 
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3 Pricing 



We price derivatives on St by taking expectation under P^.. For instance, we determine 
the price of a European option with payoff function ip : M+ — )- M and maturity t > by 

[exp i-Rt) ^(St)] . 

Special cases include: 

• Government bonds: = 1 

• Corporate bonds: ip{x) = l{a;>o} 

• Call options: ip{x) = (x — K)'^ for > 

• Power payoffs: (p{x) = x^l^^j-yoy for p e M. 

General power payoffs are not traded. But they can be priced efficiently and are helpful 
in pricing other payoffs. Moreover, for p = 1, the price of the power payoff equals the stock 
price and having a general formula for E^; [exp{—Rt)St] will allow us to obtain conditions 
for the discounted stock price to be a martingale. The case p = corresponds to a corporate 
bond. We just consider zero coupon bonds and assume there is no recovery in the case of 
default. A corporate bond with a fixed recovery can be seen as a linear combination of a 
government and a corporate bond with no recovery. So all of the arguments that follow 
can easily be extended to this case. 

In Subsection 13.11 we show how discounted moments of St can be obtained by solving 
coupled systems of generalized Riccati equations. In Subsection 13.21 we extend Fourier 
pricing methods from Carr and Madan [2] and Lee flFl to our setup. In Subsection 13.31 
we discuss the approximation of general European options with government bonds, power 
payoffs and European calls. 



3.1 Discounted moments and generalized Riccati equations 

For t>0 and x e D, define 

Ut,. := G C : E, [exp{-Rt)Sf''^'h{t<r}] < oo} 
and the discounted moments function 

ht^x{z) := E^ [ex];){-Rt)Stl{t<r}] , z G Ut^x- 

Note that Ut^x is equal to C, a half-plane or a vertical strip. Since Rt > 0, it always 
contains the imaginary axis iM.. Clearly, ht^x is finite on Ut^x and analytic in the interior of 
Ut,x- Indeed, ht^x{z) can be written as 

ht^xiz) = [exp {zst + {z - l)Rt + zKt) l{t<r}] • 

So for z in the interior of Ut^x, one can differentiate inside of the expectation to obtain 

h'^^^{z) = E^ [{st + Rt + At) exp {zst + {z - l)Rt + zAt) l{t<r}] ■ 



5 



Let 

^:R+xC"xCxC^CU {00} , 5 : IR+ x C" x C x C ^ (C U {00})" 
be solutions to the system of generalized Riccati equations 

dtA{t, u, V, w) = Go{B{t, u, V, w),v, w), A{0, u, v, w) = 0, 
< dtBx{t,u,v,w) = G{B{t,u,v,'w),v,w), Bx{0,u,v,w) = ux, (3.1) 
^Bj{t,u,v,w) = (aiY>{0jji)uj 

where the functions Go : C" x C x C ^ C and G : C" x C x C ^ C™ are given by 

Gq{u,v,w) = {au,u) + {b,u) + dv + c{w-l)+ I (e<"'^> - 1 - (uj, xj(0)) ^^K) 

Jd\{o} ^ ^ 

n 

Gi{u,v,w) = {ajU, -») + f^kjUk + SjV + 'yjjw - 1) 

fc=i 

Jd\{o} ^ ' 

The quadratic and exponential terms in Gj can cause the solution (A, B) of ( 13. ID to explode 
in finite time. However, since in condition (vii) we assumed the jumps of the process X 
to have all exponential moments, we obtain from Lemma 5.3 of Duffie et al. [6] that the 
functions Gj, i = 0, . . . , m, are analytic on C"+^. In particular, they are locally Lipschitz- 
continuous, and it follows from the Picard-Lindelof Theorem that for every (u, i), w) € 
C" X C X C there exists t* > such that equation ( 13. ID has a unique finite solution {A, B) 
for t G [0,t*). We set A and all components of B equal to 00 after the first explosion time. 
Note that the components Bj{t, ze, z - 1, z), j e J, are finite and analytic in z e C for all 
t > 0, and if Bi{t, ze, z — 1, z) is finite/analytic in z for all i e I, then so is A{t, ze,z — l,z). 
So for fixed t > 0, we define 

Vt := {z e C : Bi{t, ze, z — l,z) is finite for all i G Z}. 

and 

lt,x{z) '■= G'yiV{ze + A{t,ze,z — l,z) + {B{t,ze,z — \,z),x)} , zeVt. (3.2) 
In a first step we show the following 

Lemma 3.1. For every t > 0, Vt contains iM and hf^xi^) = k,x{z) for all x e D and z e iR. 

Proof. Fix Fx, Xn+i,Xn+2 G IR+ and denote x = (x, Xn+i,Xn+2)- Consider the process 

^ ^ j {Xt,xi + Rt,X2 + At) {ort<T 
* \ A fort>T ' 

where A is a cemetery state outside of D := x M*^"™ x to which Y jumps at time 
r. y is a Markov process with values in Z) U {A}. Since X has the properties (H1)-(H2), 
Y fulfills the assumptions of Proposition III. 2. 4 of Revuz and Yor [19] and therefore is a 
Feller process. Moreover, Rt = rudu and = Jg Xudu are of finite variation and the 
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random variable E is independent of X. So one deduces from (HI) that for all x & D and 

jE^[f{Xt,Xn+l + Rt,Xn+2 + At)l{t<T} " f (x , Xn+1, Xn+2)] 
1 

= -E^ [/(Xt,X„+i,Xn+2) - f{x, Xn+1, Xn+2)] 

+ [/(Xt, a;„+i + Rt, Xn+2 + At)e-^* - f{Xt 



^/(x) for t i 0, 



where 

Mix) 



{m + {ax,ki,xx)) 

k,l=l 



d^fix) 



+ ( 6 + /3x, V/(x) - (c + (7, xx))f{x) 



D\{0} 



dxkdxi 

(fix + (e, 0, 0)) - fix) - {S7jfix),xjim 



I JD\{0} 



ifix + (^,0,0)) - fix) - {VJu{^}fix),XJU{^}iC))) X^fiM), 



and 



• 6 = (6, d, c) 

• /3 



(n+2)x{n+2) 



/ \ 

l^jx f3jj 
6 

V 7 0/ 

By Lemma 31.7 of Sato ||20]| . the infinitesimal generator of 1" is well-defined and equal to 
A on Clib) (if /(A) is understood to be for / e CUb)). But since A is the infinitesimal 
generator of a regular affine process with values in l) U {A}, it follows from Theorem 2.7 
ofDuffie etal. lH that 

[exp ((u, Xt) + vRt + wKt) l{t<T}] = exp n, t^, w) + u, w),x)) . (3.3) 

for all it,u,v,w) G M+ x C™ x x C2_, where C_ denotes the set of all z G C with 

Re(z) < 0. In particular, for it,u,v,w) G M+ x C™ x (iM)" x C^, Ait,u,v,w) and all 
components of Bit,u,v,w) are finite. In fact, according to Theorem 2.7 ofDuffie et al. 
[6], the function B should have n + 2 components and G should have m + 2. But due to 
the special form of A, components m + 1 and m + 2 of G vanish, and the corresponding 
components of B stay equal to the initial values v and w for all times t > 0. In addition, 
Rq = Aq = 0, and one obtains ( 13. 3D . where the functions A and B solve ( 13. ID . 
Now notice that 

ht,xiz) = exp(ze)Ea; [exp (z {e,Xt) + iz - l)Rt + zAt) l{t<r}] • 

So, for all t > 0, Vi contains iM. and 

ht,xiz) = exp {ze + zs, z — 1, z) + 2: — 1, -z), 2;)} 

for all X G D and 2; G iM. □ 
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Lemma 3.2. Fix to > and zq g 14,,. Then there exists an open neighborhood W of {Iq, zq) 
in M+ X C such that 

A{t, ze, z — 1, z) and B{t, ze, z — 1, z) 
are finite for {t, z) £W as well as analytic in t and z. In particular, the set 

{{t,z) G M+ X C : z G ^t} 
is open and for all x e D, lt,x{z) is analytic in t and z. 

Proof. By Lemma 5.3 of Duffie et al. (611, the functions Gi, i = 0, . . . ,m, are analytic 
on C"+^. So, the lemma is a consequence of Theorem 1.1 in Ilyashenko and Yakovenko 
IflOll. □ 



In the following theorem we extend the identity ht^x{z) = lt,x{z) to z outside of the 
imaginary axis iR. Similar results in different settings have been given by Filipovic and 
Mayerhofer [8J, Keller-Ressel [11], Spreij and Veerman [18]. 

Denote by It the largest interval around contained in Vi n M and by the connected 
component of Vt containing 0. It follows from Lemma 13.21 that It is open in M, and it is 
clear that It ciVt- Moreover, one has 

Theorem 3.3. For all {t, x) G M+ x D, Ut^x is an open subset of C containing the strip 

{z £ C : Re{z) G It} and ht^x{z) = lt,x{z) for each z G Ut^x n V^. 

Proof. For fixed t > and x G D, one can write 

ht,x{iy) = ^x [exp{-Rt)l{t<r}] Eq, [exp{iy[st + Rt + At])] , y G M, 
where Qx is the probability measure given by 

dQx _ ew{-Rt)l{t<T} 



dFx Ex [exp(-i?t)l 



{t<r}\ 



So, up to a constant, y i-)- ht^x{w) is the characteristic function of + Rt + Kt with respect 
to Qi.. From Lemmas 13.11 and 13.21 we know that ht^x{iy) = k,x{iy) for all y G M and Vf 
is an open neighborhood of in C on which k^x is analytic. So it follows from Theorem 
7.1.1 of Lukacs I1I61I that Ut^x contains the strip {z G C : Re{z) G It} and is open. Since 
ht,x{iy) = k,xiiy) for all y G M, is connected and both functions are analytic, one also 
has ht,x{z) = lt,x{z) for z G Ut-x r\V{^. □ 

Remark 3.4. The price of a corporate zero coupon bond with no recovery in the case of 
default is given by 

Pt,, = E,. [e-^'l|5,>0}] = V'(0)- 

The price of a government bond is equal to the price of a corporate bond in a model where 
(c, 7) = (that is, the probability of default is zero). 

Corollary 3.5. The condition 

Gi(e,0,l)=0, i = 0,...,m, and ^jj = ^, (3.4) 

is sufficient for the discounted stock price exp(sf + Kt)\s^t<T} to be a martingale with respect 
to all ¥x, X G D. If all components of ej are different from 0, then ( 13. 4D is also necessary. 
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Proof. It follows from Lemma 13.21 that there exists a to > such that p G for all 
G [0, to] X [0, 1]. So if (IM1> holds, one obtains from Theorem [m and (O* that 

ht,x{l) = eMe + A{t,e,0,l) + {B{t,e, 0,l),x)) = /io,x(l) for all(t,a;) G [0,to] x D. 

Now the martingale property of exp{st + At)l{t<T} with respect to all P^, x e D, follows by 
decomposing a given interval [0, t] into finitely many intervals of length smaller than to 
and taking iterative conditional expectations. If all components of ej- are different from 
and condition ( 13. 4D is violated, there exist t < to and x e D such that /it,x(l) / ^o,a;(l)- So 
the martingale property cannot hold under F^. □ 



3.2 Pricing via Fourier transform 

In this subsection we show how to extend Fourier pricing methods from Carr and Madan 
im and Lee IflSlI to our setting. Since is used as pricing measure, one would typically 
work with models in which the discounted stock price exp{st + At)l{t<T} is a martingale 
under P^-. However, the following results just involve St for fixed t > and technically do 
not need the discounted stock price to be a martingale. 
Consider a call option with log strike k and price 



Ct,x{k) = ^OL 

Theorem 3.6. Let (t, x) G M+ x D and p > such that p + 1 G Ut^x- Then 

-pk j- -pk poo / . , \ 

ctAk) = ^ / e-'y'^gciy)dy = / Re (e-^y^gdy)) dy, (3.5) 

where 

gdy) - 



p2 _|_ p _ y2 _|_ iy(^2p + 1) 



Proof. Since Ut,x is open, there exists rj > such that p + 1 + r/ G Ut^x, and it follows from 
Corollary 2.2 in Lee IH that 

ct^k) = 0(e-(P+^)'=) for k^oo. 

In particular, 

£ (eP''ct,x{k)] 
It follows that the Fourier transform 



2 

dk < oo. 



9c{y) = / e'y''eP''ct,xik)dk 



is square-integrable in y, and one can transform back to obtain 

e^'ctAk) = ^ I e-'y^g,{y)dy. 
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This shows the first equahty of ( I3.5D . Since ct^x{k) is real-valued, one has gd—y) = gdu), 
which implies the second inequality of ( I3.5D . To conclude the proof, set z = p + iy and note 
that 



gdy) 



gSt+Ai _ ^k-Rt 



{st+Rt+Kt>k, Kt} 



dk 



st+Rt+At 



1 



z{z + 

1 



z{z + l) 



^E^. [exp{-Rt + (z + 1) log St)l{t<r}] 



□ 



To calculate prices of options with short maturities or extreme out-of-the money strikes 
Carr and Madan [2] suggest an alternative method which does not suffer from high oscil- 
lations. To adapt it to our setup, define 



dt,xik) ■■= < 



Ex 



e~^' {St - e'') 



iik < log 
if A: > log So 



Then the following holds: 

Theorem 3.7. Let {t, x) e M+ x D and p > such that 1 — p, 1 + p € Ut^x- Then 



dt,x{k) 



1 1 



where 



and 



9d{y) 



sinh(pA;) 27r 

f{y - ip) - f{y + ip) 



(3.6) 



f{y) = ^^Piii±Mfo) /.,.(0) - '^^^^hUl) + '"^ 



l + iy ' ' ■ zy ■ ■ y -ly 

Proof. Since Ut^x is open, there exists r? > such that l—p — r] and l+p + r] belong to Ut^x- 
By Corollary 2.2 of Lee UJI, 

dt^xik) = 0(e-(^'+^)l*=l) for k ^ ±oo. 

In particular, dt^x{k) and sYah.(j)k)dt^x{k) are both square-integrable in k. One easily checks 
that 

e'y''dt,x{k)dk = f{y). 



So, 



/ 



e'y^ smh{pk)dt,x{k)dk 



e^y^'"' d,x{k)dk = f(y-^p)-fiy + ^p), 

10 



and hence, 



1 1 

smh.{pk) 2tt 



-iyk 



9d{y)dy. 



□ 



There exist several extensions of the Fourier pricing methods of Carr and Madan f^. 
For example, Lee [15] shows that ( I3.5D can be adjusted in order to allow for p < and de- 
rives pricing formulas for other European derivatives, such as binary or asset-or-nothing 
options. By adjusting the proof of Theorem 5.1 of [15] to our setup, one obtains the follow- 
ing: Denote 



at,x{k) = [e-^'Stl{iogSt>k}] and bt,x{k) = [e'^'l^i^ 



■gSt>k}\ 



Then 



ct,x{k) = at,xik) - e^ht^x{k), 

and one has 

Theorem 3.8. Let {t, x) G M+ x D and p, q such that p+ I, q & Ut^x- Then 

Re (e-*^'5a(y)) dy 

qk poo 







^-pk 


at,x{k) 


= Raip) + 


vr 








bt,x{k) 


= Rb{q) + 


vr 








ct,x{k) 


= Rciv) + 


vr 



Re [e-'y'gbiy)) dy 
Re (e-'^'gciy)) dy, 



where 

9a{y) 

and 



Kx{p + ^ + w) .^ Kx{q + iy) I . 

9b{y) = — , gc{y) 



P + iy 



Raip) = < 



q + iy 

ifp<0 
Ml ifp = o, Rbiq) = { 
ifp>0 



ht,x{p+l + iy) 

p2 _|_ p _ y2 _|_ iy(^2p + 1) 



ht,x(.0) ifq<0 
ifq = 
ifq>Q 



'ht,x{l) -e^ht,x{0) ifp<-l 



Rc{p) = { 



2 





ifp = -1 
if -1 <p <0 
ifp = 
ifp > 
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3.3 Approximation of general payoffs 
3.3.1 Idea 

The prices of European options with payoff functions in the set 

Lt,x M+ M Borel-measurable such that E^. [e~'^'|99(5't)|] < 00} 

can be approximated by portfohos consisting of securities that can either be priced directly 
or with Fourier methods. For most purposes it is sufficient to use a mix of (^(0) government 
bonds, power payoffs and call options. For instance, let pi < • • • < p/ be a set of powers 
and Ki < ■ ■ ■ < Kg finitely many strike prices. Then fix s* > and determine weights 
vi,. . . ,Vf and wi, . . . ,Wghy weighted L^-regression: 



ar 



gmin / - (/5(0) - y'wisPa|^>o} - ^^^^^(s - p{s)ds, (3.7) 

\ i=i i=i J 



where p is a heuristic density approximation of St. The positive function p is meant to put 
additional weight on regions where St is more likely to lie (usually in the vicinity of [St] 
and if default is possible, around 0). If one does not have a good idea of the distribution of 
St, one can also use non-weighted regression (p = 1). 

If the integral is discretized, the optimization problem ( I3.7D becomes a finite-dimensional 
-regression. To improve the numerical stability, one can first apply Gram-Schmidt or- 
thogonalization to the basis functions 9^(0), sP'1{s>o} and {s — 

3.3.2 Example: Truncated log payoff 

We illustrate this method by approximating the price of a truncated log payoff ip{s) = 
log(s) \/ k, k e R. Note that since St = with positive probability, the truncation from 
below is crucial. 

Assume S'o = 1 and k = —1. We consider three ways of approximating ip with linear 
combinations of 101 instruments: 

1. A government bond and power payoffs of powers 0.05, 0.1, . . . , 5 

2. A government bond and call options with strikes 0.03, 0.06, . . . , 3 

3. A government bond, power payoffs of powers 0.1, 0.2, ... ,5 and call options with 
strikes 0.06, . . . , 3 



We let s* = 3. As heuristic density for St we use p given by 

p{x) 



exp(— lOx) X < 0.5 

exp(-lOlx-ll) 0.5<x<1.5. (3.8) 
_exp(— 5) X > 1.5 



This choice of p assigns more weight to and points around 1. But it is just an example of 
what one could choose. Depending on the model and the value of t one may want to use 
different functions p. 
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0.5 
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1 1.5 2 

Stock price 

Approximation with government bonds, power payoffs and calls 



o 



-0.02 




1 1.5 

Stock price 



Figure 1: Comparison of different approximation methods. Note that the scales are dif- 
ferent! 



Figure [T] shows the errors of the three approximation methods. It can be seen that for 
the truncated log payoff methods 2 and 3 give a much better approximation than method 
1. The errors of methods 2 and 3 are similar. But since prices of power payoffs are easier 
to calculate than those of call options, method 3 is significantly faster. 

Log payoffs are useful in the pricing and hedging of variance swaps on futures. Let 
Fu, < u < t,he the price of a futures contract on St and consider a variance swap with 
time-t cash-flow 

where = to < ■ ■ ■ < tj = t are the discrete monitoring points (usually daily), K is the 
strike and C is a cap on the payoff (typically C = 2.5K). If F„ is modeled as a positive 
continuous martingale of the form dFu = a^FudWu, then the probability of hitting the cap 
C is negligible and one can approximate the sum with an integral: 

~ 7 / ^ld.u - K. 
t Jo 

It has been noticed by Dupire (H and Neuberger IfTTll that 

^ aldu = 2 (^^ j^dFu - log St + logFo^ , 
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and therefore, 



-E, [log Fo- log St] 



In a diffusion model with the possibihty of default the cap is crucial since it is hit in 
case of default. If one neglects the probability that the cap is hit before default or that 
log St < k := log(Fo) — t{C + K) /2 if there is no default, one can approximate T,t as follows: 



{r>t} 



aldu -K\+ l{r<t}C 



tAT Y \ 2 

dFu - {log{St) V A;) + logFo - K - l{r<t}T 



F, 



^u- 



-dFu. 



In the special case where the default intensity is constant, the expectation of the last 
integral is zero, and one can price according to 



E, [Sj] « -E, [log Fo - log{St V k)\ - K. 



4 Hedging 

In this section we consider a subclass of affine models in which European contingent 
claims can perfectly be hedged by dynamically trading in sufficiently many liquid se- 
curities. Assume that condition ( 13. 4D is fulfilled. So that the discounted stock price 
exp(sf + Af)l|i<^} is a martingale under allPi.,x G D. It is well-known that it is impossible 
to replicate contingents claims with finitely many hedging instruments in a model where 
the underlying has jumps of infinitely many different sizes. Therefore, we here require 
the jump measures v and Hi, i G Z, to be of the following form: 

M 

u = ^ Vq5y^ for f q > and different points yi, . . . , yAf in G ^\{0} (4.1) 

q=l 

Mi 

in = ^ Vig6y^^ for Vig > aud different points yn, . . . , yiMi in G D\{0}, (4.2) 

q=l 

where 5y^ , 6y^^ are Dirac measures and for M = or Mj = 0, or /xj are understood to be 
zero, respectively. The integrability condition (vii) is then trivially satisfied. 

In the following theorem we are going to show that the process X has a realization as 
the unique strong solution of an SDE of the from 

dXt = biXt)dt + a{Xt)dWt+ [ k{Xt-,z)N{dt,dz), (4.3) 

where W is an n-dimensional Brownian motion and an independent Poisson random 
measure on with Lebesgue measure as intensity measure, a : D W^^"^ has to be a 
measurable function satisfying 

aa^{x) = a + ^ XiUi for all x e D, (4.4) 
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and the functions 6 : D — ^ M" and k : D x D are of the following form: 

b{x) = b + (3x- / XjiO^idO - yZ^i / Xju{i}{Ol^i{dO^ 
Jd\{0} ^ Jd\{0} 

{Vq ifVq-l<Z<Vq 
yiq if Vm + Xj^jM-i + XiVi(^q_i) < Z < Vm + E}=1 ^j^jM, + XiViq , 

ifz>VM + ET=i^J^jM, 

where Vq = Ylp=i^p Viq = X]p=i ^jp- Theorem 14. II is an extension of Theorem 8.1 in 
Filipovic and Mayerhofer |83. Its proof is given in the appendix. 

Theorem 4.1. If v and ^i, i ^Z, are of the form ( I4.1D -( |4^ . then there exists a measurable 
function a : D —> M"^" satisfying ( 14. 4D such that the SDE 04. 3D has for all initial conditions 
X e D a unique strong solution. It is Feller process satisfying (H1)-(H2) corresponding to 
the parameters (a, a, b, /?, m, /i). 

In general, one needs L + 1 := n + M + YlT^i + 2 instruments to hedge all sources 
of risk. We assume the first one to be a money market account yielding an instantaneous 
return of rt. In addition, one needs one instrument for each of the n components of the 
Brownian motion W, one for each of the possible jumps of X and one for the jump to 
default (without default, L instruments are generally enough). Of course, if the hedging 
instruments are redundant, not all European contingent claims can be replicated; precise 
conditions for a given option to be replicable are given in ( 14. 5D below. The hedging in- 
struments should also be liquidly traded. For instance, in addition to the money market 
account, one could use a mix of instruments of the following types: 

• Stock shares 

• Government bonds 

• Corporate bonds 

• Vanilla options. 

All of these can be viewed as European options with different payoff functions. Let us 
denote the set of hedging instruments different from the money market account by 

$ = {(tl,V9i),...,(tL,(^L)}, 

where ti, . . . , are the maturities and Lpi e Lti^x the payoff functions. 

Now consider a European option with maturity t < min {ti, . . . , t^} and payoff function 
(f G Lt^x- At time 0, its price is 

7r(t,x) := [exp(-i?t)(/?(S't)] , 

and for u G [0, t], 

_ J Tr{t-u,Xu) ifn<r 

"^""l E4e-(«*-^") I X„] (^(0) ifn>r ' 

After the default time r, Cu behaves like a government bond and can be hedged accord- 
ingly. To design the hedge before time r, we introduce the following sensitivity parame- 
ters: 
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• Classical Greeks: for g = 1, . . . , n: 

d 

• Sensitivities to the jumps corresponding to v. for all g = 1, . . . , M: 

Jt,x = '^{t,x + yq) -Tr{t,x) 

• Sensitivities to the jumps corresponding to ^: for all i G Z and q = 1, . . . , M^: 

• Sensitivity to default: 

Dt,^ = [exp(-i?t)^(0)] - [exp(-i?t)(^(St)] . 

Example 4.2. Consider a power payoff (p{s) = for some p > such that G Lt^x- Then 
the sensitivity parameters are given by 

d 

= ht,x (p) = Bg{t,p£,p-l,p) ht,x {p) 

Jt!x = ht,x+yM - ht,x{p) 

Dt,x = -ht,x{p)- 

Example 4.3. For a European call option with log strike k and maturity t, the classical 
Greeks are 

Hlx = Q^^^Ak) = -;r I (e"'"^-.5c(y)) dy, 

where 



dx„9c{y) = d: 



1 p2 j_ p _ y2 iy{2p + 1) 

^ Bq{t, {p + l + iy)e,p + iy,p + l + iy)ht,x{p + l + iy) 
p'^ + p-y^ + iy{2p + 1) 

for some p> such that 1 + p G Ut,x- The other sensitivities are 

Jt,x = Ct,x+y,ik) - Ct,xik) 
Jt% = Ct,x+yig{k) - Ct,x{k) 

Dt,x = -ct,x{k). 
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As shown in Examples 14.21 and 14.31 the sensitivity parameters of power payoffs and 
vanilla options can be given in closed or almost closed form. The payoff 93 e Lt^x can be 
approximated with a linear combination of government bonds, power payoffs and Euro- 
pean calls of maturity t as in Subsection 13. 3[ The sensitivities of ip are then approximated 
by the the same linear combination of the sensitivities of the government bond, power 
payoffs and European calls. If X is a solution of the SDE ( 14. 3D and vr is a C^'^-function on 
(0, t] X D, one obtains from Ito's formula that for u <t At, 



dCu = "drift" du+ Hl^,x^_agg,{Xu-)dW^' + Dt-u,x^^dl{,^^} 

q,q' = l 

+ / {TT{t-u,Xu- + k{Xy,-,z)) -7r{t-u,Xu-))dN{du,dz), 
Jr+ 

which can be written as 

n 

dCu = "drift" du+ Hl^,x^_agg>{Xu-)dW^' + Dt-u,x^_dl{r<u} 

q,q' = l 

M m Mi 

+ E J?-u,x.-dm + E E Ji^u,x.-d^:^, 

q=l i=l q=l 

where N'' and N'''^ are Poisson processes with stochastic intensity depending on X but 
independent of the payoff function ip. For all / = 1, . . . , L, define vr' and C' analogously to vr 
and C, and assume the functions vr' are all C^'^ on (0, t] x D. To hedge before default, one 
has to invest in the hedging instruments such that the resulting portfolio has the same 
sensitivities. That is, one tries to find i9(n, x) € such that for all < u < t A r. 



'Hl^,x. = Eti^\u,X^)Hll\^ 



Xit 



TTtl 

^t-u,X^ 



n 

'^t-U,Xu 



'^t-u,X^ 



Ei=i^\u,X^)jl;\ 



T.i=i^\u.X^)j\fXx^ 



(4.5) 



T.t=i^\u,x^)j[t-u,x^ 



'^t-u,x — 1^1=1^ \^^^u)'J^_u^Xu 

where on the left side are the sensitivities of ip{St) and on the right, indexed by I, the sensi- 
tivities of the hedging instruments (f\Sti). Xu) is the number of hedging instrument 
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I in the hedging portfoho before default while the amount Cu - J2iLi '&\u, Xu)Cl is held 
in the money market account. Since e^^^'Cu and e~^"C^, / = 1, . . . , L are all martingales 
under P^, this is a self-financing strategy replicating C until time t A r. If default happens 
before time t, one simply holds (f{0) zero coupon government bonds from then until time 
t. 

( I4.5D is a system of L linear equations with L unknowns. It may not have a solution if 
it is degenerate. But if the family $ of hedging instruments is such that ( 14. 5D has full rank 
for all < li < t A r, then any European contingent claim can be replicated by dynamic 
trading. 



5 Heston model with stochastic interest rates and possibil- 
ity of default 

As an example we discuss a Heston-type stochastic volatility model with stochastic inter- 
est rates and possibility of default. It extends the model of Carr and Schoutens lUl and 
can easily be extended further to include more risk factors. Let {Xt)t>Q be a process with 
values in D = M.\ x M moving according to 



dX^ = Ki {6i - X}) dt + i]iJxldWl (5.1) 



= K2{e2-X'^)dt + r]2^X^dW^ (5.2) 
dXf = -]-Xldt+ JxldWf (5.3) 



2 

for a 3-dimensional Brownian motion W with correlation matrix 




and non-negative constants ki, K2, 9i, 62, r/i, r/2. Since X^ and X"^ are autonomous square- 
root processes, the system ( I5.1D -( |5^ has for all initial conditions x G x R a unique 
strong solution, and it follows as in the proof of Theorem 14.11 that it is a Feller process 
satisfying (H1)-(H2) with parameters 




p?7i \ / 

• a = 0, = i I , = i 

1/ V 



-Ki 

• 6 = (ki01,K2^2,O), /3= I -K2 

-1/2 

Note that one cannot have correlation between W"" and W"^ or W"^ and without de- 
stroying the affine structure of X. So we introduce dependence between s, the volatility 
process X^, r and A by setting 

St = xl n = d + 61x1 + 82x1 \t = c + ^1x1 + ^2x1 

Rt = Tudu, At= Xudu 
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for non-negative constants c, 71, 72, d, 61,82. Then 



St = exp {st + Rt + At) l{t<T} 



satisfies the SDE 



dSt = St-(^[c + d+ (71 + Si)Xl + (72 + ,52) X^] dt + ^dWi - c?l|,<t}^ 



5.1 Pricing equations and moment explosions 



It follows from Corollary 13.51 that the discounted price exp (sj + At) l{t<T} is a martingale. 
The discounted moments function is of the form 



ht,x{z) 



exp(^(t, (0, 0, z), z-l,z) + {B{t, (0, 0, z),z - 1, z),x)) 
exp ( A{t, z) + Bi{t, z)xi + -B2 (t , z)x2 + 2x3 ) , 



where 



dtA(t, z) = Ki0iBi{t, z) + K292B2{t, z) + (c + d){z - 1) 

dtBi{t, z) = \r]lBl{t, z) + {priiz - Ki)Bi{t, z) + {^z + 71 + 6i){z - 1) 

dtB2{t, Z) = \liBl{t, Z) - K2B2{t, z) + (72 + 62){z - 1) 

i(0,x) = ^1(0, z) = ^2(0,2) =0. 



(5.4) 



In this special case, Bi and B2 are both solutions of scalar Riccati ODEs that can be 
obtained explicitly. The explosion times of the discounted moments 

t*{p) := sup {t>0: E[eM-Rt)Sfl{t<r}] < 00} 

can also be determined in closed form: 

(arctan ^ + 7rl{^<o}) if n < 

ilogH±v^ iiu>0,v>0, 

cxD otherwise 

where 



t*{p) 



(5.5) 



PViP - 1^1, 
^P + 7i + 6i] {p- 1). 



The derivation of 05. 5D is analogous to the derivation of the same formula for the Heston 
model which can be found, for example, in Andersen and Piterbarg [1]. Figure [2] shows 
the decay of t*{p) for p ±00. A sample plot of the implied volatility surfaces for the 
two cases of positive probability of default and no default is given in Figure [H For the 
calculation of implied volatilities we used the yield on a government bond as the interest 
rate in the Black-Scholes formula. As one would expect, credit risk contributes towards a 
higher implied volatility, especially at longer maturities or at extreme strikes. This effect 
can help explain why implied volatilities usually exceed realized volatilities. 
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0.34 



0.32^ 




Figure 3: Implied volatility surfaces with positive probability of default (upper surface) 
vs. no default (lower surface). Parameters: ki = 0.06, K2 = 0.04, i]i = 0.2, r]2 = 0.1, 9i = 1, 
6»2 = 0.3, p = -0.6, c = 0.02 (no default: c = 0), 71 = 0.01 (no default: 71 = 0), 72 = 0.01 (no 
default: 72 = 0), d = 0.01, 61 = 0.1, 62 = 0.1, = 0.05, = 0.03, X^ = 0. 
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5.2 Hedging 

Denote by A'^{t,z), Bl{t,z) and B2{t,z) the solutions of the Riccati equations ( I5.4D for 
c = 7i = 72 = (no default). Then the price of a zero coupon government bond with 
maturity t is given by 

= exp(iO(t, 0) + ^0(t, 0)xi + 0)X2). 

Assume there exists at > such that 

d^,Cu,Ak)B^iu,0) - d,,CuAk)BUu,0) + (5.6) 
[B0(n,0)B2(^x,0) - ^0(n,0)^i(7x,0)][a,3C„,,(A:) - c„,,(A:)] 

for all < n < t and x ^ D. Then every European contingent claim can be hedged if the 
following four instruments can liquidly be traded: the stock, a zero coupon government 
bond, a zero coupon corporate bond and a call option with log strike k, the latter three all 
with maturity t. 

Indeed, by Remarks 14.21 and l4.3[ the hedging parameters are as follows: 
. Stock: = H^^^ = 0, H^i = D^, = -e^'^. 

• Government bond: H^:^ = Bl{u,Q)Pl^, H^;^ = B^{u,0)P^^^, H^:^ = 0, = 0. 

• Corporate bond: H^'x = Bi{u,0)Pu,x, H^Jx = B2{u,0)Pu,x, H^Jx = 0, 1?^^ = -Pu,x. 

• Call option: Hu;x = d^^Cu,x{k), g = 1, 2, 3, D'^^^ = -Cu,x{k), 
and ( 15.61 1 is equivalent to 

e'^Pl,Pu,x{-dx,Cu,x{k)Bl{u,Q) + d^,Cu,x{k)B\{u,Q) 
+[^0(u,0)^2(n,0) - B0(tx,0)^i(n,0)][5,3C„,,(A;) - c„,,(A;)]} / 0. 

So the matrix 

/ ^0(^z,0)pO, Bi{u,{))Pu,x d^,cu,x{k) \ 
Bl{u,Q)Pl, B2iu,0)Pu,x d,,cu,x{k) 
e"^ dxsCu,x{k) 

\ — e^^ ~Pu,x ~Cu,x{k) / 

has full rank, and the system ( 14. 5D always has a solution. 



A Proof of Theorem 4.1 



It is shown in Theorem 8.1 of Filipovic and Mayerhofer fSl that there exists a measurable 
function a : D ^ M"^ " satisfying (gS such that the SDE 



dXt = b{Xt)dt + a{Xt)dWt 



(A.1) 
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has for each initial condition x & D a unique strong solution X^^\ Now set tq := and 
define iteratively 

Tq := inf|t>Tg_i: [ [ k{xi''_^^\ z)N{du,dz) o] (A.2) 
I Jo Jr+ J 

where y*^"?^ is the solution of the SDE dA.ll ) on [r^, cxo) with initial condition 

Y(^) = X(f 1) + k{xilZ'\z)Nirq,dz). 

Since X^'^^^^ is RCLL and the intensity measure of is Lebesgue measure, it can be seen 
from ( IA.2I ) that Tq > Tq^i a.s. So the process 

is the unique strong solution of ( I4.3D on [0,roo), where Tq. It remains to 

show that Too = oo and X^"^^ is a Feller process satisfying (H1)-(H2) with parameters 
(a, a, 6, /3, m, ^u). To do that we introduce the counting process 



By Ito's formula, 



tAT„ 



Jo 

is for all f G C''^{D x M+) and g > 1, a martingale, where 

fc,/=i ^ ' fc=i \ /=i / 



+ / {fix + z + 1) - /(x, z) - {Vjfix, z),xjm) HdO 
Jd\{o} 

m „ 

+ 5^ / + + -/(a:^,^) - (VJ7u{^}/(a;,^),X:7u{^}(0))a;^M^(^^0• 

~:T Jd\\o\ 



lD\{0} 

But ^ is the infinitesimal generator of a regular affine process X with values in D x M+. So 
by Theorem 2.7 of Duffie et al. (Bl, C^{D x R+) is a core of A, and it follows from Theorem 
4.4.1 of Ethier and Kurtz Il7l that the martingale problem for A [c|(DxIR+) is well-posed. 
Moreover, the stopping times Tq are exit times: 

Tq = mf{t>0: Zt^ [0,(7-1/2)}. 

Therefore, one obtains from Theorem 4.6.1 of Ethier and Kurtz |[7ll that the stopped mar- 
tingale problem corresponding to A |c|(DxR+) and D x [0,g — 1/2) is well-posed for all 
q. Hence, (X^^-* , Z(/^^^)f>o has the same distribution as {Xt/\f^)t>o, where fq is the g-th 
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jump time of X. Since X is RCLL and Vm + yiMi^if^ is the jump intensity of Z, 

we conclude that almost surely, the process xj:"^^^ jumps at most finitely many times on 
compact time intervals. In particular, Tqc = oo, and Theorem 14.11 follows from Theorem 
2.7 of Duffie et al. [6| since the first n components of X form a regular affine processes 
with infinitesimal generator A acting on / € C^{D) like 



kA=l 



^/(^) = J2 ("'^^ + i^^M^^x)) + {b + (3x, V/(x)) 



{fix + - /(^) - (V^/(x), 

D\{0} 
m „ 

+ E / (/(^ + 0- fix) - {VJu{^}f{x),XJU{^}{0)) X,lli{di). 
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